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Abstract 

We show how to represent a class of expressions involving discrete sums over partitions 
as matrix models. We apply this technique to the partition functions of 2* theories, i.e. 
Seiberg-Witten theories with the massive hypermultiplet in the adjoint representation. 
We consider theories in four, five and six dimensions, and obtain new matrix models 
respectively of rational, trigonometric and elliptic type. The matrix models for five- and 
six-dimensional U(l) theories are derived from the topological vertex construction related 
to curves of genus one and two. 
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1 Introduction 



Matrix models have become an important ingredient of theoretical physics over the last 
two decades. The information about the matrix model is encoded in the partition function 



where the measure DM represents an integral over an appropriate set of N x N matrices 
M. The solution of the matrix model is encoded in the spectral curve which arises from 
eigenvalues condensing in the extrema of the potential V in the large N limit. Taking at 
the same time h — > with the 't Hooft coupling g s N = const, the solution of the matrix 
integral is given by coefficients F g which arise in perturbative expansion in h of F = log Z. 

In this paper we show how to represent a class of expressions involving sums over 
two-dimensional partitions as matrix models. Such expressions often arise as partition 
functions of Af = 2 supersymmetric gauge theories and topological strings. In particular 
we will apply this technique to so-called 2* theories, i.e. Seiberg-Witten theories with the 
massive hypermultiplet in the adjoint representation. 

To start with we recall some other relations established between matrix models and 
supersymmetric field theories, and explain how our approach differs from them. First of 
all, one such relation is the Dijkgraaf-Vafa correspondence [TJ. According to this corre- 
spondence the low-energy effective superpotential of Af = 1 theory is computed by the 
matrix model, whose potential is identified with the tree level superpotential. While im- 
mense amount of work has been done to verify and apply this correspondence in the Af = 1 
context, it is also possible to extract from it some information about Af = 2 theories. To 
do that one has to consider Af = 2 theory broken to Af = 1 by a superpotential treated as 
a perturbation, adjust the superpotential in a way which allows to choose a relevant point 
in the Coulomb branch, apply the Dijkgraaf-Vafa conjecture for the Af = 1 situation and 
eventually identify certain parameters and set the deformation to zero to get Af = 2 data 
[2]. In this way the prepotential Tq of Af = 2 theories and its gravitational correction T\ 
have been obtained from matrix model calculations in some cases [HI [HE]. 

The results of Dijkgraaf of Vafa allow us to get the information related to F-terms 
in Af = 1 theories due to underlying holomorphicity. Another sector of these theories is 
encoded in D-terms whose full determination would require going beyond the holomorphic 
regime. Nonetheless just the holomorphic information is sufficient to solve more rigid, but 
still non-trivial Af = 2 theories [HI Ej - Seiberg-Witten solution of these theories is encoded 




(1) 
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in a family of Seiberg-Witten curves. The research in this direction culminated in the 
derivation of the Nekrasov partition function 0, E] for Af = 2 theories. These partition 
functions Z^ekrasov arise from the localization and are given by sums over two-dimensional 
partitions which label instanton configurations. They can be expanded as 

Z^ekrasov C^-'s—Q 3 

and apart from the prepotential fF Q they also encode the whole series of gravitational 
corrections T g . 

Gauge theories with Af = 2 supersymmetry are also intimately related to string theory 
on local Calabi-Yau manifolds via geometric engineering [10|. In particular partition func- 
tions of the large radius A-model topological strings on toric Calabi-Yau manifolds agree 
with Nekrasov partition functions for 5-dimensional gauge theories upon appropriate iden- 
tification of parameters, and reduce to 4-dimensional result in a proper limit [TT 1 [T2 l [TB I [T4] . 
As the problem of computing the large radius topological string partition functions on toric 
manifolds is essentially solved by the topological vertex theory [TH US], so are the corre- 
sponding gauge theories which can be geometrically engineered. In this correspondence 
the Seiberg-Witten curve is encoded in the target space of the mirror B-model geometry. 
More generally, each toric threefold encoded in the diagram constructed from topological 
vertices is dual to a geometry based on a corresponding Riemann surface. 

One immediately notices similarities between matrix models, Af = 2 theories, and 
topological strings. In each system there is a Riemann surface which governs the solution: 
the spectral curve, the Seiberg-Witten, and the curve in the B-model target space. It is 
advantageous to relate explicitly gauge theories and topological strings to matrix models, 
both conceptually, as well as with a hope to apply techniques from one side of such an 
equivalence to the other side. The Dijkgraaf-Vafa correspondence can be regarded as one 
such possibility. However this allows to obtain only JF and T\ data from the corresponding 
matrix model F and Fx, while for g > 2 gauge theory JF 9 's do not agree with matrix model 
F s 's [HE]. Moreover this lowest genus Af = 2 answer is not explicitly given by the matrix 
integral in this case, but it arises only after certain identifications of parameters. We also 
note that the idea of c-deformation was advanced pEE] to relate the all genus partition 
function of the matrix model, including non-planar diagrams, to the superpotential of the 
deformed gauge theory, however this is not immediately helpful from the matrix model 
point of view on gravitational corrections to the Af = 2 theories. 

Another connection between matrix models and gauge theories or topological strings 
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is provided by the formalism of Eynard and Orantin. In [19] they showed how to solve 
the loop equations and derived a recursive procedure to get correlation functions and the 
partition function of a matrix model. In fact, this recursion and its solution is determined 
by the spectral curve and the associated differential alone. One can therefore start with 
any curve and a differential and derive the corresponding set of symplectic invariants, 
without the knowledge of the matrix model, or even without the underlying matrix model. 
It was indeed checked in [20] that such invariants, computed for a Seiberg-Witten curve or 
a B-model curve, agree with gauge theory and topological string invariants. However this 
relation to matrix models is not explicit, and the lack of the matrix model potential still 
prevents one from considering phenomena such as eigenvalue tunelling. 

In this paper we take yet another approach and present how partition functions of a 
certain class of gauge theories and topological strings can be written as matrix models. 
A first step in this direction was taken in [21], where a matrix model for bundles over P 1 
was written down. This result was confirmed in |22| , where in more general context the 
Plancherel measure was rewritten as a matrix model. In [23] these results were generalized 
to pure SU(n), Af = 2 theories in 4 and 5 dimensions. In 5-dimensional case it was also 
possible to take into account Chern-Simons terms, which gave rise to more complicated 
toric geometries and multi-matrix models. 

The matrix models which we derive here and those in [23] in principle encode all 
gauge theory amplitudes fF g for arbitrary g, without a need of taking any special limit or 
identification of parameters. The price for that is quite involved form of the matrix model 
potentials: while to get the correct Tq and T\ in the Dijkgraaf-Vafa case it is sufficient to 
choose a polynomial of order (n + 1) (to make eigenvalues condense in n extrema and give 
rise to the genus (n — 1) Seiberg-Witten curve), the potentials which we derive abound in 
polylogarithms and their quantum generalizations. Of course these new potentials have the 
correct number of extrema to give rise to the curve of the appropriate genus (it was checked 
explicitly in [23] that the matrix model spectral curve agrees with the Seiberg-Witten curve 
for the pure SU{2) theory). 

In a wider context, in this paper we present a more general method of turning into 
matrix models certain expressions given by sums over partitions. This method is based 
on the (last line of the) formula ((4)) . We use this method in a particular case of Af = 2* 
theories. The matrix models which we derive for such theories in various dimensions are 
given by formulae (fPOl). flT3l. (l20l). and fl25l). 

The plan of the paper is as follows. In section [2] a construction of matrix models for 
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a class of expressions involving discrete sums over partitions is presented. In section [3] 
partition functions of 4-dimensional 2* theories, obtained from instanton sums, are refor- 
mulated as matrix models. In section 0] a matrix model is derived for a partition function of 
the 5-dimensional theory obtained from the topological vertex calculus and related to the 
genus one curve. In section [5] we derive a matrix model for 6-dimensional theory, related 
to the genus two curve. Section [6] contains conclusions and discussion. 



2 The hook length formula and matrix models 

In this section we explain how to construct matrix models encoding discrete sums over 
partitions, with arbitrary dependence on the hook lengths. While the general idea of 
turning discrete sums into matrix integrals is the same as in [22l I23|, here we generalize it 
to a wide class of expressions. Such expressions arise in particular as partition functions 
of supersymmetric theories and topological strings. In the following sections we apply this 
method to the case of Af = 2* theories in various dimensions. 

To start with we fix some notation. For a partition R = (Ri, R2, • • • , Rn), with Ri > 
R 2 > . . . Rn > 0, we define 

fH = Ri - i + N, (2) 

which constitute a strictly decreasing sequence. The hook length of a box □ = G R 
(with i and j denoting respectively the row and the column of a box) is defined as 

h(D)=R i + R t j -i-j + l, (3) 

where R* is a transposed matrix. 

We now present a formula which is a basis of many calculations in what follows. For a 
function cp such that <p(n) ^ for n / we have 



n l tt <p(Rj - Rj - i + j) rr jfjhj - hj) 

(p(h(D)) -1-1 <p(j-i) -1-1 <p(j-i) 

□ei? rv v " i<i<j«x> ' \<i<j<oo rw ' 

n ^fe¥)(nn 1 



fU - i) yv fJi fe = i v(. k -j + N). 

N hj 

( n *-y)(nn^)' w 

j=l k=l^ K ' 



l<i<j<N 
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The equalities in the first two lines were derived and used (among the others) in [T2l[T3]. 
This is however the third line which immediately allows to relate expressions involving 
products over hook lengths with matrix models, and which we use repeatedly in next 
sections. To go from the second to the third line we notice that the first product in the 
denominator can be rewritten as Yli<i<j<N < P(j ~~ *) = n^i* ' ' '^(j)]- Then each 
expression in the square bracket combines with one of the products \\ k 3 = \ <p{k — j + N) 
into the product of the form Ylk=i I n total we get N such expressions, each one for 

fixed j = 1,...,JV. 

We are interested in deriving matrix models corresponding to the sums over partitions 
involving products of the above type. To find such models we first truncate sums over all 
partitions to sum over partitions of a fixed number of N rows. Then we use the formula 
(BJ to replace the dependence on hook lengths by the dependence on hi defined in (j2J). The 
first factor in the third line Y[ ^p(hi-hj) turns into some generalization of the Vandermonde 
determinant and becomes a part of the integration measure, while the product Ylk=i ~^{kj 
gives rise to the matrix model potential V (after writing it as an exponential). The product 
T[f=i represents the trace. We rewrite the above sum as continuous integrals by introducing 
an auxiliary function with poles in all integer values of hi and treating Ui = hhi (or 
equivalently Ui = —g s hi in five- and six-dimensional cases) as continuous matrix eigenvalus 
in the 't Hooft limit 

h -> 0, N -> oo, hN = const. (5) 

Because of the rescaling of hi, the dependence on h (or g s ) appears not only in front of 
the potential term in (00) but it also enters the potential in a nontrivial way. We therefore 
expand the potentials in a series in h (or g s ), in particular in order to identify the lowest 
order terms which survive the 't Hooft limit. 

Example 1: for R — (6, 3, 1, 1) one can check the formula (jlj) explicitly, substituting in 
this case (hi, hi, h 3 , h±) = (9, 5, 2, 1) and the following hook lengths of all boxes: 
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Example 2: for ip(x) = x 2 the product over hook lengths reproduces the Plancherel 
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measure [22]. The equality Q implies then 

Dei? v ' ' h u ...,h N l<i<j<N j=l r 

Introducing Ui = hhi allows to write this as a (continuous) matrix model, with the Van- 
dermonde determinant arising from the product Y\(hi — hj) 2 , and the potential given by 
the logarithm of the T function (which arises as the analytic continuation of the factorial) . 
This way we immediately recover the results of [22] . 



3 4-dimensional theories 

In this section we derive a matrix model for 4-dimensional 2* theory. This theory can 
be constructed from M = 4 theory by adding mass m to the hypermultiplet. This theory 
has a microscopic coupling r , and the associated parameter Q = e 2mr ° can be used in 
instanton counting. We also use the following notation: 

Q = e~ T , m = hfx, Ui = hhi. (7) 

In general in what follows hi will be defined as in (j2J), while Ui will become a continuous 
parameter after taking the limit (J5J). 

First we consider U(l) theory, or more precisely its twisted version, for which the 
partition function localizes on instanton configurations labeled by partitions R. It was 
shown in [9] that this partition function can be written as a sum over partitions 

R U&R y ' 

In this case each summand is a product over all boxes of R and the dependence on each 
box is only via the hook length. Identifying <f^\x) = x 2i 2 we can now use our hook 
length formula. The first factor in Q becomes 

n f u u\ h t - h i) = n n-SAr^ =fi hi ~ h ' 



,< t<1 <N i<t<i<N < ,! - - h > - A (hi-hj+p) y h, - hi - f, 

Using n£ =1 (fc + A 4 ) = T( r(p+i) ,> we ^ n< ^ ^ ne contribution to the potential 



h 



n 



k 2 - /i 2 T(l + h + n)T(l + h- n) 

^ ~ r(i + / i)r(i-^)r(i + /^- 
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To replace sums by integrals over matrix eigenvalus, we introduce the function [22l [23] 

f(x) = -xT(-x)T(x)e-™ : 



siniirx) ' 



which has simple poles at all integer values of the argument. To find expansion of the po- 
tential in h we also make use of the expansion of the digamma function (i.e. the logarithmic 
derivative of the T function) 

= i log r(*) = -i + log x - jr (9) 

m=l 

where B?. m are Bernoulli numbers. 

Combining the above, and writing the second factor in equation (@J as an exponent, 
we obtain the matrix model of the form (pQ) with 



N N 

Uj — U 



vmt = n^n^i <">> 



Ui — Uj-m 

k=l i^j J 

Tr u(i), m u2 u — m 

Ud \ u ) = Tu + u log — + m log h 

u — m z u + m 



Jl\ u2 y B2kh2k ( 1 1 — 

2 ° g u 2 - m 2 ^ 2ife(l - 2k) \(u + m) 2k - 1 (u - m) 2k - 1 u 2 ^ 1 



As expected the matrix model measure turns out to be a mass deformation of the 
Vandermonde determinant. The terms in the first line of the potential are the only terms 
essential in the planar limit and are shown in figure Q3 The terms in the last line are 
proportional to higher powers of K and vanish in the limit (jSJ) ■ 

Now we consider a generalization of the above result to U(n) theory. A particular 
point on the Coulomb branch we are interested in is specified by the values ai = hpi and 
we introduce = ai — a&. The instanton part of the partition function for U(n) theory 
with the adjoint hypermultiplet reads [HI [9] 

R=(RM,...,R<~ n )) 

z i-T a ik + h(Rj l) ~ Rf ] aik + m + h(j - i) 

n (mMm aik + ^' - l) <*» + ™ + KRf ] - Rf 

This partition function involves a summation over n partitions R^ and explicitly is 
not of the form (111). To deal with it we perform the same trick as in [23]. For each partition 
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V(u) 




u 



Figure 1: Planar limit of the matrix model potential for 4-dimensional, £/(!), 2* theory. 



R( k ) we introduce h\ defined as a slight generalization of ((2j), 

hf ] = R l -i + N + Pk , 

and then concatenate 

U _ f h W /.CO 1,(2) /,( 2 ) 1>) />)x hn\ 

ftj=i,...,nAf — l^i ,...,n N ,n x ,...,n N , , % ,...,n N ). [lz) 

This effectively leads to the one-matrix model with eigenvalues Ui = hhi. After manipula- 
tions similar to those in [23], we obtain the matrix model of the form (Dp) with 



nN N 

Uj — U 



vmt = n^n^s us) 



Ui - Uj - m 

k=l i^j J 



^ « = Tu+ y\(u- ai) log- - + mlog ■ + 

L (u — aiY — m z u — ai + m 



oo 



2HT 



' 2&(1 — 2fc) V (m — a; + m) 2fc 1 (u — a/ — m) 2fc 1 (w — a/) 



2&-1 



ft (w - a;) 2 

+ 77 lOE 



2 (it — a^) 2 — m 2 . 



Again the measure is a mass deformation of the Vandermonde determinant, and only 
the terms in the first line of V^ n \u), shown in figure [2, contribute in the planar limit. 

In principle one could generalize the above 4-dimensional matrix model to 5 and 6 di- 
mensions by replacing the instanton factor in f fTTl ) by its trigonometric and elliptic counter- 
parts [H1CE2]. This would lead to the corresponding generalizations of ([TBI . In what follows 
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we will consider such generalizations only for U(l) theory, using the results obtained from 
the topological vertex computations and stressing the relation of these theories to genus 
one and two curves. 

V(u) 




Figure 2: Planar limit of the matrix model potential for 4-dimensional, U(3),2* theory. 



4 5-dimensional U(l) theory 

In this section we consider 5-dimensional U(l) theory compactified on a circle of size (3, 
taking the advantage of the topological vertex computations presented in [12]. Following 
the topological string convention, we use the parameter g s instead of h, as well as 

q = e~ 9 % q h > = e - 9sh \ Ui = -g s h h 

with Ui referring to the continuum matrix model counterpart of hi defined in (j2J). We also 
use the following notation for the deformed quantities 

h 

[h] = q- h ' 2 - q h '\ [h]\ = l[(q- l/2 - q l/2 ). 

i=i 

The instanton part of the partition function is related to the (p, q) five-brane web or 
the topological vertex diagram shown in figure [3l Two Kahler parameter of this diagram, 
T and T m , relate to the gauge coupling r and the hypermultiplet mass m, so that 

Q m = e~ Tm = e^ m , Q = e ~ T = Q~ 1 e 2mT0 . (14) 
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Figure 3: Geometric realization of the 5-dimensional, 2*, U(l) theory in terms 
of the periodic topological vertex diagram. Slashes on horizontal legs denote 
their identification along a circle. 



The instanton partition function encoded in the diagram in figure [3] is given in terms 
of the topological vertex as 

P,R 

with the amplitude Cpqr derived in [15]. This partition function was manipulated in p2] 
to the form 

7 ud) _ sr nw TT ^-Qmq h(U) )ii-Q m q' h{U) ) h 

z 5d 11 n _ ah(P))(i _ a -h(n)\ 

R DeR K H JK * ' 

This expression is explicitly of the form (j3J) . If we define 

{ X ' Q) ~ (1 - Qq*)(l - Qq-»y {lb) 
a straightforward computation gives 

TV • i Ui—Uj 

<P S } ) (h i -h j ,Q n ) = e > II -7 Ul -l-T m ■ (17) 

l<i<j<N i^j blllLl 2 

This gives rise to the measure of the matrix model we are after. 

To derive a matrix model potential we use similar tools as in [23]. We write the 
quantum dilogarithm as 

oo 

g(x) = Y[(l-x- 1 (f). 



i=l 
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It vanishes g(q h ) = for h a positive integer, and at such points its derivative is 

(-, \2 mh n -h(h-l)/2 

q'(q h ) = - 9{ ' Q . 

9{q ) q h (l-q h ) 9 (q- h ) 

Therefore the following function has simple poles with residue 1 for x = q h with h 6 N 



(7(1; e 9= 6 e 2 ss 



(1 — x)y/xg(x)g(x~ 1 ) 
Let us also note that for x = q l 

-1(1+1) /4 gOO 



g(x~ 1 q p ) 2 1 

exp(— log(xg" p )log(xg" p+1 )). 



([Z-p]!) 2 ^(1)2 ^2]ogg 

Also note [22] that to get the expansion of the potential in g s one needs 



^ OO ~Q y 

\ogg{x) = V Li 2 - m {x~ 1 )^f 

q s i ml 



and in consequence, for q h = e u , we have 

JJ(1 _ Q q k ) = *® = e'Ts ^=o (Li 2 - m (Q)-Li 2 _(^)) (lg) 
fc=i S'Iq^J 
Combining the above, the potential for 5-dimensional theory arises from the following 
contribution to the formula (TJj) 

TT 1 = n h 9{Qm)9{Qm) -j- Em=0 (2Li 2 - m (e")-Li2- m (Q m e")-Li 2 - m (Q- 1 e")) 

(19) 

To conclude, rewriting the expression (fT5l) using the above manipulations, and up to 
overall constants, we get the matrix model expression of the form (pQ), with K replaced by 
g s , and with 



N N ■ i Ui—Uj 

= n^n TtTT. < 2 °) 

jfc=l i^j blllU 2 

= -u(T + T m ) + 2U 2 (e u )-U 2 (Q m e u )-U 2 (Q^e u ) + 
+g s u + % log 



2 & (l-Q m e«)(g-Q-ie«) 

00 75 2fc 

+ E-(^-( 2Li 2(^)( eM ) -Li2(i- fc) (Q m e«) -Li 2(1 _ fc) (Q-V) 



12 



As expected the measure of this matrix model is the trigonometric and massive defor- 
mation of the Vandermonde determinant. The first line of the potential V^ l \u) contains 
the only terms which contribute in the planar limit, which are shown in figure 0J 

V(u) 
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Figure 4: Planar limit of the matrix model potential for 5-dimensional, £7(1), 2* theory. 



5 6-dimensional U{1) theory 

Generalization of the U(l), 2* theory to 6 dimensions can be understood in terms of 
the toric diagram with two periodic directions, representing the curve of genus 2, as shown 
in figure El In addition to parameters Q and Q m introduced in f fT4l) and relevant for the 
situation in figure [3l now there is an additional Kahler parameter 

Qi = e~ T \ 

The rest of the notation is the same as in section 01 

Instanton partition function for 6-dimensional theory with U(l) gauge group, related 
to the above genus two surface, is given by the topological vertex amplitude 

zT = X)(-Q) |JI| (-W P| (-Qi) |5| ^(?)^^(g)- 

P,R,S 
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Figure 5: Geometric realization of the 6-dimensional, 2*, U(l) theory in terms 
of the doubly periodic topological vertex diagram. Slashes on horizontal and 
vertical legs denote their identification along a circle. 



This expression was manipulated in |12| to the following form 

7 um Vr)Wn (i-Q m g fe(D) )(i-Q m g- ft(D) ) x m , 

z 6d ~ 11 n _„/»(□))(! -n-m) [ ' 

r neR v H A H ' 

■q (i - Q k P Q m g h{n) )( i - Q k p Q m q- h(D) )(i - Q k p Q^ l q h{D) )(i - QjQmV MD) ) 



;i-Q*gMO))2(l-Q* g -A(n))2 



fc=l 
where 

Qp QlQm- 

This again involves summands of the form Q which allows to make contact with matrix 
models. We define 

^q*){l-Q k p Q^q- 

(22) 



(7(1)/ x_ (7(1)/ n \ "1 T l x ) \± ^ipH ) 

<P Gd W - <P U ^ W 11 (1 _ QkpQmqX){1 _ Qk Qmq - x){1 _ Q*Q- l5 -)(l - Q*Q-} q -*) ' 



with Q m ) given in ((THJ). 

The measure of the 6-dimensional matrix model arises from the familiar by now prod- 
ucts over y^j (hi — hj), which after some manipulations become 

11 ^ n 3 ) e 11 m< t*-h s \Q ) [I6) 

l<i<j<N i^j L ^ mH \^P> 

with the notation 

e 1 (t\Q) = i J2(-i) n Q^ + ^. 
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The manipulations necessary to get the above measure include in particular replacing 
infinite products over k in by the theta function using the following form of the 

Jacobi triple product identity 

0i(t\Q) 



IK 1 - <A)(i - Q kt ^ 1 ~ Q k ) = -Qi /8(t i /2 _ ,-i /2) 

The matrix model potential arises from the terms of the form 



where we introduced the quantum polylogarithm 

W*.Q)= E fl -OTini« - (24) 



Special cases of this function include the quantum dilogarithm L 2 considered in [24], as 
well as the quantum trilogarithm L 3 which appears in (25j EH [27] . 

Collecting the above results, and after some more algebraic manipulations, we obtain 
the 6-dimensional matrix model 



= n^n^^ (25) 



fc=i t^j 

oo 



^V) = -«(T + T m )+^ + ^^(2Li 2 _ fc (e u )-Li 2 _ fe (Q m e")-Li 2 _ fe (Q- 1 e u ) 

fc=0 ^ 



-2L 3 _ fc (e u ,g p ) +L 3 _ fc (Q m e u ,Q p ) +L 3 _ fc (Q- 1 e u ,Q 



pi 



The measure of this potential is the expected elliptic and massive generalization char- 
acteristic for compactified 6-dimensional theories [321 [27]. In the planar limit only the 
terms in V§ d (u) proportional to g® ~ 1 contribute, i.e. only the linear term and the first 
term for k = from the infinite sum, which involves the dilogarithms Li 2 and quantum 
trilogarithms L 3 . 



6 Discussion 

In this paper we derived matrix models for 2* theories in 4, 5 and 6 dimensions starting 
from expressions for partition functions obtained from instanton counting or topological 
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vertex calculations. The measure of these matrix models is respectively of rational, trigono- 
metric and elliptic type and includes also the massive deformation. The potentials are quite 
involved, and their expansion in h (or g s ) involves ordinary and quantum polylogarithms. 

Our matrix models can also be regarded as complementing the results of |T2J . In 
that paper matrix models, topological vertex calculations and elliptic genera were related 
to each other, also in the context of 2* theories. The matrix models discussed there are 
however of Dijkgraaf-Vafa type and are obtained starting from 1* theory, i.e. Af = 1 theory 
with three adjoint chiral superfields $1=1,2,3 of mass m and the tree-level superpotential 
Tr($i[$ 2 ) $3])- As explained in the introduction in this case there are difficulties in relating 
higher genus matrix model F g to gauge theory amplitudes T a . In contrary, by the very 
construction, the matrix models for 2* derived here explicitly encode all gauge theory T g . 

It is natural to expect that the method we used can still be generalized. An interesting 
feature of the matrix models presented here is the fact that the formula (j3J) implies that 
both the measure and the potential of the matrix model are encoded and derived from the 
same function ip. To obtain matrix models in yet more general situations would require 
relaxing this feature. One aim of this program would be to convert into matrix models 
topological string partition functions for local Calabi-Yau manifolds encoded in arbitrary 
toric diagrams. 

On the other hand, the use of the formula Q is not limited to the gauge theory or 
topological string partition functions, and with its help matrix models could be derived 
for other interesting physical and mathematical systems. Some such applications related 
to the simplest case of the Plancherel measure in the formula (j§J) were discussed in [22] . 
As for the next-to-simplest case of the 4-dimensional 2* theory (jHJ) let us just note, that 
some interest aroused in the field of enumerative combinatorics PSJ [29] due to the formula 
found by Nekrasov and Okounkov in [9] 

00 

zT = IK 1 - W 1 - 

n=l 

It would be interesting to see how combinatorial features of this formula relate to the 
combinatorics of diagrams of the matrix model ( TTOl) . 

Our results should also be of interest from several other points of view. On one hand, 
Seiberg-Witten theories with matter, including SU(2),2* theory in 4 dimensions, were 
considered recently in [30] from the viewpoint of holomorphicity and modularity. Also the 
relation to the formalism of Eynard-Orantin was discussed there. Relating matrix models 
derived here to those results would be an interesting check of all these approaches. 
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On the other hand, Af = 2 theories in various dimensions are related to integrable 
systems. In 4, 5 and 6 dimensions theories these are respectively Calogero-Moser system 
[3T] . Ruijsenaars-Schneider system [32], [33] , and their generalization considered in [T21 15~4~|. 
Recently the new viewpoint on the relation to integrable systems was given in [27]. It 
would be of interest to understand relations between matrix models derived here and 
those integrable systems. 

2* theories with U(l) gauge group in 5 and 6 dimensions, which correspond to toric 
geometries representing curves of genus one and two, were also considered recently in the 
context of the refined vertex [35] and related to cylindric partitions [36, 37J. One might 
wonder if there is any way of obtaining such a refined information from more general matrix 
models. 

Finally, the genus two curve relevant for 6-dimensional U(l) theory arises also in the 
context of dyon counting in J\f = 4 string theory [38l 15*5] . The dyon counting function 
is given by the Igusa cusp form of weight 10, and in [30] it was shown to arise also as 
the genus one contribution to the partition function for the genus two curve. Finding any 
relations between both these systems, as well as between matrix models derived here and 
©-modules analyzed in [""I"]"""""], is certainly worth more effort. 
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